(N 
O 
O 
(N 

> 
O 

oo 
(N 



(N 



Oh! 



X 



IKDA-01/22 
SI-2001-7 
hep-ph/Ollom? 
October 2001 



Deep Inelastic Scattering 
with Tagged Initial State Radiation: 



(N , 

; Complete 0{a) leptonic QED Corrections 



Harald Anlauf|J| 

^ . Institut fiir Kernphysik, Darmstadt University of Technology, 

O^'. 64289 Darmstadt, Germany 



and 

Fachbereich Physik, Siegen University, 57068 Siegen, Germany 

Abstract 

In this paper we extend the calculation of the QED corrections to 
deep inelastic lepton-proton scattering with a tagged photon, taking 
into account the full corrections on the lepton side. Comparing to 
previous results that were obtained by considering only large logarith- 
mic terms at leading and next-to-leading accuray, we find that the 
difference is in general quite small, however, it may be significant in 
the region of large y and small x. 
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1 Introduction 



One of the major aims of the experiments at the HERA ep colhder is the mea- 
surement of the structure functions of the proton, F2(x, Q^) and Fl^XjQ"^), 
over a broad range of the kinematic variables. Especially the domain of small 
Bjorken x < 10^^ and momentum transfer of the order of a few GeV^ 
and below is of particular interest, as it provides a challenge for attempts 
towards a complete, quantitative understanding of the dynamics of quarks 
and gluons inside the nucleon. 

Of these structure functions, the longitudinal one, Fl, is much more diffi- 
cult to access. There exist several ways to separately extract F2 and Fl from 
the experimental data. The most obvious, direct method requires to run the 
collider at different (i.e., lower) center-of-mass energies, which may not be 
desirable from the point of view of other parts of the physics program. 

Indirect methods usually require substantial input from theory and de- 
pend more or less on modeling the hadronic final state, like extrapolations 
or QCD fits (see e.g., or the measurement of the azimuthal angle 

distribution of final state hadrons, as suggested by Gehrmann However, 
these methods can be used at fixed collision energy. 

Another direct method was suggested by Krasny et al. and utilizes ra- 
diative events with an exclusive hard photon registered in the forward photon 
detector (PD). Such a device is actually part of the luminosity monitoring 
system of the HI and ZEUS experiments, and will continue to exist after the 
HERA luminosity upgrade. The idea of this method is that emission of pho- 
tons in a direction close to the incoming electron corresponds to a reduction 
in the effective beam energy. The effective electron energy for each radiative 
event is determined from the energy of the hard photon observed (tagged) in 
the PD. 

Besides measuring F^, radiative events extend the accessible kinematic 
range to lower values of Q^. The potential of this method is supported by 
preliminary results from the HI collaboration of an analysis at low for F2 
|^,P,^,|T^ (for earlier analyses that did not take into account QED radiative 



corrections see |lTT|,p!2|). The feasibility of the corresponding determination 
of Fl was studied in However, with currently analyzed data sets it is 
not yet possible to compete with F^ from extrapolations or QCD fits . 

A precise analysis of experimental data requires the inclusion of radiative 
corrections. The most important ones are the QED corrections on the lepton 
side that have been discussed at the leading logarithmic level and 
taking into account next-to- leading logarithms As the difference 

between leading and next-to-leading logarithmic terms could be quite signif- 
icant, reaching of the order of 5 percent in some regions, it appears difficult 
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to estimate the remaining uncertainty due to non-logarithmic terms. It is 
the purpose of the present work to close this particular gap. 

We may focus on the QED corrections on the lepton side since they 
form a gauge invariant subset of the full corrections to the process under 
consideration. Furthermore, one expects that the QED corrections on the 
hadronic side are significantly smaller, as the typical hadronic mass scale of 
the order of the proton mass is much larger than the electron mass, so there 
is no comparably large logarithm. In fact, emission of photons off quarks can 
essentially be absorbed into suitably defined parton distribution functions. 



and the net effect turns out to be numerically negligible . 

The outline of the present paper is as follows. After introducing our 
notation in section 2, we devote section 3 to a discussion of the individual 
contributions to the leptonic corrections with forward photon tagging, and 
apply these corrections in section 4 to the case of radiative deep inelastic 
scattering. In section 5 we present some numerical results and our conclu- 
sions. The appendices collect several formulae that are useful for a numerical 
implementation of our calculation. 



2 Born cross section 

A a starting point, let us introduce our notation in the context of the lowest 
order contribution to the radiative, semi-inclusive deep inelastic scattering 
process 

e{p) + p{P) e{p') + 7(A:) + X(P') , (1) 

where the photon 7 is assumed to be measured in the photon detector. 

We choose invariant kinematic variables that use the measured, scattered 
electron and take into account the energy loss due to photon emission 



X 



2P-{p-p' - k) 



P-ijp-p' -k) 
y = P-ip-k) • 

For HERA conditions, the polar angle ■(9^ of the tagged photon (as mea- 
sured with respect to the incident electron beam) will be very small, -^^ < "f^o, 
with -(^o being about 0.45 mrad in the case of the PD of HI. We shall also 
assume below that the scattering angle of the electron, 6, is always much 
larger than -^q. Therefore, the energy fraction of the electron after initial 
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state radiation reads 



2P-{p- k) Ee - 



S 



E. 



xyS 



(3) 



where E^ is the electron beam energy, E^ represents the energy deposited in 
the forward PD, and S = 2P ■ p. 

The relation between the shifted variables and the standard Bjorken 
variables of deep inelastic scattering reads: 



Q 



2 Q' 



X 



xy 



y = 1 - z{l ~ y) 



(4) 



z 1- z{l-y) ' 

The Born cross section, integrated over the solid angle of the photon 



detector (0 < "i?^ < ''^O; m/ E^ <^ 1) takes a factorized form. 



1 dVeo 



a 



with 



y dxdydz 27r 



P(z,Lo)S(x,y,Q2) 



Q'^xy'^ 



2(l-y)-2xV 



M2 



+ il + Ax' 



y 



2 



(5) 



(6) 



and 



1 + z' 



2z 



In 



a 



1 - n(-Q2 



R 



R{x, 



l + 4x^^ 



^M^\ E2{x,Q'' 



1 



(7) 



2xFi(x,g2; 

The obvious advantage of the use of shifted variables (H) is their direct 
appearance in the argument of the structure functions in the Born cross 
section. 



3 Leptonic corrections 

In the expression for the lowest order cross section, eq. (^), we encountered 
the logarithm Lq =: In^o of the quantity 

Co = ^ . (8) 
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Although i^o <^ 1, we have Co ^ 1 and Lq ^ 6.5 ^ 1 for the conditions of 
the HERA PD's. We shall therefore consistently neglect contributions that 
are of order (9('i9o) or (9(Co"^)- 

As explained in the introduction, the subset of leptonic QED corrections 
is gauge invariant, and it also factorizes, thus allowing a discussion isolated 
from the hadronic part. Besides keeping things more transparent, this also 
facilitates reusing the results in other calculations. 

In the present section, we shall therefore consider the Compton subprocess 

e(pi)+7*(-g) ^e(p2)+7(fc) , (9) 

with the emission angle of the photon being integrated over the PD, while 
taking into account the corrections from virtual and real QED corrections. 
While performing this integration, we require that the remaining part of the 
amplitude for the full process (i.e., 7* + p — > X) depends only weakly on the 
small transverse momentum of the forward photon. 



3.1 Compton tensor 

Let be the matrix element of the Compton scattering process (P), with 
the index /i describing the polarization state of the virtual photon. Adopting 
the notation of |2^, we define the Compton tensor 

= Mf-^'\M:^'-^''r . (10) 

spins 

Using current conservation, this tensor is conveniently decomposed as 
follows: 

Pi^y = 9f.u {Bg + + ^ Pi^Pj^ {B,j + ^Ti^ , (11) 

_ Qf^Qu - _ _ Pi ■ q ■ _ 1 9 

The expressions for the quantities Bij corresponding to the Born approx- 
imation are: 

Bg = ^4(s + M)2 + (t + M)2] -2mV(^ + ^) , 



Bu = — — , B22 = — — , B12 = B21 = , (12) 



8m^ Aq^ Sm? 

~ ^ ! B22 = — ^ 

St S2 St t2 

s = 2^2 ■ k , t = —2pi ■ k , 
u = (pi-p2)^, s + t + M = g^. 
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Note that for almost coUinear emission, ~ (1 — z)pi, we may neglect 
the transverse momentum of the emitted photon in the tensor decomposition 
( ]TT| ) and use momentum conservation to set p2 = zpi. 

The kinematic variables of the Compton subprocess are related to those 
of the radiative DIS process via: 

u = -—, = {pi-k-p2f ^ -Q"^ , s^- -Q"^ . 

z z 

The quantities T^, in (0) denote the radiative corrections to the 
Compton tensor. 



3.2 Virtual and soft corrections 



The virtual corrections to the Compton tensor, as calculated in [^], are 
conveniently decomposed into a piece containing the universal infrared sin- 
gular contributions, which are proportional to the Born contributions, and 
an infrared finite remainder: 



(13) 



where 



A 



7r2 9 



p = 41n-(LQ - 1) - + 3Lq + 31n2; + - - - , Lq = In ^ . (14) 



m 



The parameter A in the above expression is a fictitious photon mass regulating 
the IR divergency. 

Performing the integration over photon angles, we obtain 



{-QhiJ^v + ^zpi^^piy) X 
1 



(15) 



1 - z 



i + ^p)p{z.l,)-^t\+o{^Iq') , 



2tx' 



where 



T 
A 

B 
C 



{A\nz + B)P{z,Lo) + CLo + D 
2LQ-Lo-21n(l-2) , 

ln^2 - 2Li2(l -z) - - , 



2z 

T^z 



\nz — z 
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1 _ fir + 4^2 

D = -— (U2(l-z)+lnzln(l-z)) 

1 — z 

,9/ X 82; , 47r^ 

- 2z\n%l-z)-\ \nz z + 1, 

^ ' l-z 3 



X 

Li2(x) = - ! ^ln(l-i/). 

J y 



The single and double logarithmic terms in Lq and Lq of the above expression 
agree with ref. [|17|. Details of the calculation will be given elsewhere pT[ . 



The dependence of the virtual corrections on the unphysical parameter A 
is canceled by the contribution from emission of an additional soft photon, as 
usual. Requiring that the energy fraction of the second (soft) photon in units 
of the energy of the incoming electron does not exceed e, with e <C 1, and 
adding the contribution from soft photon emission to the virtual correction 
then amounts to the replacement of the quantity p in (|TB|) by p, see p]| : 



p = 2{Lq - 1) In^ + 3Lq + 31nz - In^F - ^ - ^ + 2Li2 ] , (16) 



with 



E'e 

Y = — ^ and c = cos^ = cosZ(p, p") (17) 



being the relative energy of the scattered electron and the cosine of the 
scattering angle in the lab system, respectively. 



3.3 Double hard bremsstrahlung 

For the case of double photon emission, we define the 'double Compton ten- 
sor' analogously to (|T0|) as 

= ^ E Mf-^'^^M:^*-^'^^r , (18) 

spins 

where now M'p*^'^''^^ is the matrix element of the double Compton process 
process 

e{pi) + l*{-q) ^ e{jp2) + l{ki) + 7(^2) , (19) 

with the index p describing the polarization of the virtual photon. 

For the kinematic invariants of this subprocess we shall use the notation: 

Zi = 2pi -ki , z'- = 2p2 ■ ki , 
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a = 2ki-k2 = {h + , 

A = -[{pi - ki - k2f' - m^] = zi + Z2 - a , 

A' = [{p2 + ki + k2f -m^]= z[ + z'^ + a , 

Ql = -{Pi - P2)^ = 2{pi ■ p2 - m^) , 

Ql = -q'^ = Qf + Zi + Z2 - z[ - z'^ - a . 



3.3.1 Double collinear emission 



When both photons are emitted almost coUinearly to the incoming electron, 
[ki ~ Xipi, k2 — X2P1), we may neglect the transverse momenta of the 
photons with respect to the incoming electron, and the double Compton 
tensor takes a simple form: 



K 



2-coll 



4 [-^M^Q? + 4:Z {pi^pi^ 

1+^2 1 z 

A2 



X 



X1X2 Z\Z2 

m" I z^ 



^ - + - + 



- 2 



-Az- 



X2zf 



+ 



2 

X ]^ 



m 
A2 



Here ri 



Xi, r2 



1 1 
— + — 

Zl Z2 
1 - X2, 



m 



A2 

and z = 1 



X1X2 V ziA 

(1 - z){rir2 + z) 

X1X2Z1Z2 

1 1^^ 
— + — 

Zl Z2 



rf + zr2 ^^^2 + ZTi 



Z2A 



(20) 



Xi 



X2- This expression is 



consistent with Merenkov where leading and next-to-leading logarithms 
were calculated. 

Performing the integration over the photon angles and over the relative 
photon energies, taking into account the symmetry factor 1/2! for the emitted 
photons, one obtains: 

' e)Q{x2-e)5{xi + X2-(l-z))Kl-^'^'' (21) 



— / dA;idA;2 0(xi 



[-9t^uQl + 4:Z {pi^piu 



X 



a 



p(2) p(2),IR-div. p(2),IR-fin. 
log ' nonlog nonlog 



The previously known leading terms containing double and single logarithms 

(2) 

Lq are contained in -Piog - The nonleading terms are split into an infrared 
divergent piece that depends on Ine, and an infrared finite piece. An outline 
of the calculation and expressions are given in appendix ^. 

It is worth to mention that the terms depending on the soft-photon cutoff 
parameter e in (^1]) do factor nicely, as expected from the usual soft-photon 
factorization: 



[-~gi.uQl + ^z{pi^Pi,)\ X (^|-)'p(z,Lo)-2(Lo-l)ln^ 



:22i 
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3.3.2 Final state coUinear radiation 



Consider now the kinematic region where one photon, say, photon 1, is emit- 
ted almost coUinearly to the incoming electron, i.e., ki ~ xip, and the other 
close to the outgoing electron, so that k2 ~ ^(p2 + ^2)- The double Compton 
tensor then simplifies to 



FSR 



Xi Zi 

i + (i-0' 1 



X 



-2(1-0:1) 



m 



X 



exhibiting the expected complete factorization of coUinear initial and final 
state radiation, respectively. 

Integrating this expression over the emission angles of photon 1, one ob- 
tains 



v/ 



dQi K 



FSR 



2 

Xi 

2 



P{l-Xi,Lo) X 

2 



+ 4(1 - xi) ■ {pi^piu) 



X 



^ ^2 ^'2' . 



Note that we have not yet integrated over the angles of the final state 
photon. The treatment of final state radiation depends on the experimental 
situation, i.e., whether the detector is able to resolve a photon coUinear to 
the electron, or whether it just measures the sum of their energies. 



3.3.3 Semi-collinear emission 

The final case covers the kinematic range where one photon is emitted almost 
coUinearly to the incoming electron, while the other photon is emitted at 
an angle ■i?2 > 1^0, but not coUinear to the final electron. We denote this 
kinematic domain as the semi-collinear one. 

In order to be consistent with the above calculation of the double coUinear 
emission, we shall perform the angular integration over the coUinear photon 
and drop all contributions of the order 0{j)i^ and (9((^^^). We find indeed 
factorization of initial state radiation for large emission angles of the second 
photon, i.e., 'd^ ^ ^o- However, in the vicinity of the forward cone that is 
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defined by the solid angle of the PD, there are contributions from further 
terms with a complicated Z2-dependence that spoil a naive factorization. 
These additional terms fall off rapidly and essentially contribute only in the 
small region < ^2 ^ "^^o- 

Assuming that photons in this narrow region outside the PD will not be 
measured, we integrate these terms over angles and split their contribution 
schematically as follows: 

^ [ dn, [ dn^ K';r~"'''' 



/?2 r 

~ — / dr^a 

1?2>1?0 



~, Ar,{Qf + Z2)r + {nQ^ - z',y 



^ P{n,Lo] 



riZ2Z2 riZ2Z2\ XiTi 

+ [-dtivQf + 4(ri - X2) (pi/xPii.)] H{xi, X2) , (23) 

X\X2 

with ri = 1 — Xi- The explicit expression for the function H{xi,X2), which 
collects the mentioned non-factorizing, quasi-collinear terms, is given in ap- 
pendix ^ It is infrared-finite and does not contain any logarithm of a large 
scale. 

The integrand in the first part on the r.h.s. of ( p3l) can be rewritten as 

lp(ri,Lo)-- 

Xi Ti 



-P(r^,Lo)--B^rirm,P2), (24) 



where for the sake of consistency one should drop terms of order on the 



r.h.s. For a discussion and further details see ref. pT |. 

Since in our decomposition of phase space only photon 1 reaches the PD, 
we have to identify ri by z and xi by 1 — 2; in the above expressions. However, 
we still need to integrate over the phase space of the other photon that is 
emitted at large angles. This calculation depends on the complete scattering 
process and in general requires a numerical integration. 



4 Radiative DIS 

Let us now turn to the description of the radiative scattering process 
e(p)+p(P)^e(p')+X(P0+7(^i) ( + 7(^2)) . 
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As argued in the introduction, the leading raditiative corrections to the pro- 
cess stem from emission of photons off the lepton hne.Q It is now straightfor- 
ward to contract the radiatively corrected Compton tensor of the previous 
section with the hadron tensor 

H^,u{P,qh) = 47i(-g^,Fi{xh,Ql) + P^,Pu-J—F2{xh,Ql 

\ ■ qh 

= ATi(^-~g^,F,{xh,Ql) + P,P,'^F2{xh^Ql)^ , 

P — P -P " qh _ Qh ^2 _ 2 

— qhv 2 ' ~ TTn ' ^h ~ qh ■ 

qi 2P ■ qh 

Here qh denotes the four-momentum transfer to the hadronic system. Taking 
into account the emission of two photons, qh = Pi — P2 — ki — k2. 

Applying the results from the previous section, we find for the contribu- 
tion from virtual and soft corrections to the cross section: 

where p is taken from ([T6| ) with 

y = ^ = z{l-y)+xy^ , 

c = COS ^ = ^ . (26) 

z[ \ - y)Ee + xyEp 

In the calculation of the contributions from the emission of two hard 
photons, we decompose the phase space into three regions discussed in the 
previous section (see also fl^): i) both hard photons hit the forward photon 
detector, i.e., both are emitted within a narrow cone around the electron 
beam {{}i^2 < "^^o, "^^o ^ 1); one photon is tagged in the PD, while the other 
is collinear to the outgoing electron (-i^g = Z(A;2,p') < "^q)', and finally in) 
the second photon is emitted at large angles (i.e., outside the defined narrow 
cones) with respect to both incoming and outgoing electron momenta. For 
the sake of simplicity, we assume that m/i?e -C -C 1. 

The contribution from the kinematic region i), with both hard photons 
being tagged in the PD, but only the sum of their energies measured, reads: 



y dxdydz Svr^ 



p(2) p(2),IR-div. p(2),IR-fin. 
log ~'~ nonlog ~'~ nonlog 



(27) 



"'^ Another gauge invariant set of large corrections is due to the running of the QED 
couphng a. Note that we have aheady included it in the expression for the Born cross 
section, so we won't discuss it further. 
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see eq. (pT]). 

In region ii) we need to distinguish between the cases of whether a pho- 
ton emitted close to the outgoing electron can be detected separately (exclu- 
sively), or whether its energy and momentum are measured together with the 
electron (inclusively), as this affects the reconstructed kinematic variables. 

For the exclusive event selection, when only the scattered electron is de- 
tected, we obtain 



1 d3 



a. 



77 



^,excl 



a 



y dx dy dz Air'^ 



P(z,Lo) I 



1 + 



L - 1 



:28) 



where L = ln(E^¥[)Vm2) = In {E^^'Q^m^) + 2 In F, and E/ = S(x/, yf, Q)) 
is an implicit function of C, via the relation between the "internal" kinematic 
variables Xf,yf,Qj- and the "external" ones x,y,Q'^ (see refs. [|16|,0,[T^ for 
more details). The integration limits explicitly depend on the method for 
the determination of kinematic variables. 

In the case of a calorimetric event selection, where only the sum of the 
energies of the outgoing electron and coUinear photon is measured and taken 
into account in the determination of the kinematic variables, the correspond- 
ing contribution reads 



1 d^alZ 



;,cai a 



y dxdydz Air"^ 



P(z,L 



0) 



1 + c! 
i-c 



L- l + 21nC +1-C 



(29) 

see also [p!6|,p!7|,p!^]. 

For the contribution from the semi-coUinear region in) we apply the re- 
sults from section 13.3.31 to obtain: 



1 dV 



3^77 



a 



y dxdydz 



+ I dx2 H{1 - z,X2)T.{xt,yt,Q 



t) 5 



(30) 



47r2 J z — X2 

e 

where the 'radiation kernel' P and the boundaries of integration for the log- 



arithmic part are given in |T^. For the second, quasi-coUinear contribution, 
we have used the abbreviations 



Xt 



{z — X2)xy 
zy-X2 



yt 



zy - X2 

z -X2 



, Qt=Q 



2Z - X2 



(31) 
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and the upper limit of integration is 



X2 



zy 



(32) 



The total contribution from QED radiative corrections is finally obtained 
by adding up (|25D , (pT]), (pOf), and, depending on the chosen event selec- 
One easily verifies that the unphysical IR regularization 



tion, dl) or (|2|) 
parameter e cancels in the sum 



5 Results and Discussion 

In this section we shall present numerical results obtained from the above 
expressions and compare to next-to-leading radiative corrections |T^. As 
input we used 

Ee = 27.5 GeV , Ep = 820 GeV , ^0 = 0.5 mrad . (33) 



We chose the ALLM97 parameterization as structure function with 
-R = 0, no cuts were applied to the phase space of the second photon, and 
we assumed a calorimetric event selection. Furthermore, we took a fixed 
representative angular resolution of -^'q = 50 mrad. 
Figure |1| compares the radiative correction 

dV 

5rc = -rr, 1 , (34) 

d^Born 

calculated with next-to-leading logarithmic and full 0{a) accuracy for the 
electron method at x = 10~^ and x = 0.1 and for a tagged energy of E-pD = 



20 GeV. Similar to the well-known QED corrections to DIS (see e.g., |24 
and references cited therein), the corrections are large and positive for large 
y, while they may become large and negative for ^ — >■ at large x. 

Furthermore, one sees that the difference between the full result and the 
one at next-to-leading logarithmic accuracy is very small; it is typically 



at the permille level except for large y, where the phase space for the 'lost 
photon' gets large, see (0), and for small x. In this case it easily reaches 
values of the order of a percent. This may be important in view of the 
increased statistics expected at the upgraded HERA collider. 

Looking at the individual contributions beyond the next-to-leading log- 
arithmic approximation, one finds a significant cancellation between terms 
from virtual+soft corrections and from double coUinear emission. Their sum 
is typically of the order permille, and it depends only on z, c.f. (^) and (P?]). 
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y 

Figure 1: Radiative corrections (5rc (iH) with full leptonic 0{a) vs. next-to- 
leading logarithmic accuracy at x = 10""^ and x = 0.1 and a tagged photon 
energy of 20 GeV. No cuts have been applied to the phase space of the 
second (semi-coUinear) photon. 



13 



Therefore the total difference can be mainly attributed to the contribution 
of photons emitted into a small region outside but close to the PD. 

To summarize, we have calculated the full leptonic QED corrections to 
DIS with tagged initial state radiation. We find that the corrections beyond 
next-to-leading logarithmic accuracy are typically very small, but can still 
be significant, i.e., of the order of a percent, in the interesting region of small 
X, where the use of radiative events appears most promising. 



A Integrals for double collinear emission 



The calculation of the contribution from double collinear emission assumes 
that only the sum of the photon energies, (1 — z)E, will be measured in the 
forward photon detector. 

We split the integration of expression (^1]) over the restricted two-photon 
phase space in the following way: 



dA:i dk2 e{^o - A) &{^o - A) 6 ((1 - z) - {xi + X2)) 
1 



xi dxi / X2 dx2 5 ((1 - z) - {xi + X2)) 



(35) 



In the last line we adopted the notation of Arbuzov et al. [p5]. 



(36) 



for the angular part of the integrals. 



A.l Integrals over photon angles 



We shall now provide the relevant angular integrals needed for the contri- 
bution of two photons emitted almost coUinearly to the incoming electron. 
The calculation is performed under the assumption that E'^d'l/m? ^ 1 and 
For details we refer the reader to 
Using the abbreviation 



1^1 < 1. 



21 



Lo 



we obtain: 



^12^2 



X1X2 



^0 ; 
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4A 



Z2 



x\x2ri 
1 

—i- 

XiX2rj { 2 



Lo + ln 



X2ri 
1-z 



In 



ri(l - z) 



X\Z 



^L2 + Loln^+ Li2 
2 xxz 



X2 



lni^-l + 



+ 5 ( cosV^; 



Xxr2 
X2ri 



XiZ 



z 



+ Lis 

+ 



X2 



X\Z 

r\r2 — 22; 



+ S cos-;/'; 



a;ir2 



(37) 



z 



H2x2) + ''^''~^'^ ln 



1-z 



ri(3r2-4:z). rir2 + 4:XiZ. . 
H m 2; m ri m r2 



2z 2^ 

Az — rir2 . , 
H m (77 + X2ri + Xir2 cos ip) 



2z 



2z 



rir2 

~2r 



In {rj + Xir2 + X2ri cos ■0) 



where 2; 



1 - a;i - 0:2, ri,2 = 1 - 2:1,2, and 

, 1 2xiX2 

cos -0 = 1 , 

rir2 

77 = {xi + .T2)(xi + a;2 - 4x1X2) , 

I, f VTT2txTx^ + tx + i 



E{t;x) = |W 



+ Li2 
+ Li2 



2 

(1 + t),x 



(38) 



Vl + 2tx 



tx 



{l~t)x 



Vl + 2tx + x2 + tx + l, 
The integrals that contribute only non-enhanced terms (i.e., no I/q's) are: 





1 


_2;i2;2A_ 


^1^2 



(1 — 2;) ln(l — 2;) + 2; In 2; — ri In ri — r2 In r2 
Xi In xi — X2 In X2 



1 ^^ ri(l-^) 



X1X2 



2^2 



1'*'2 



X12; 

xiy ^^ ri{l - y) 



X2 



xiy 
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(39) 





1 







xl{3 — 1x\) In ■ 



X\Z 



+ X2(3 — 2x2) m 



2:22; 



+ In- 



rir2 



2xiX2 



The remaining integrals can be obtained from those given above by exchang- 
ing photons 1 and 2, i.e., X\ ^ X2 and ri ^ r2- 

It should be noted that the coefficients of the double and single logarith- 



mic terms (Lq and Lq) agree with refs. ||22|,P5 



A. 2 Integrals over relative photon energy 

Having performed the angular integration, we still need to integrate over the 
relativ photon energy, see (0). As this integral will be infrared-divergent, 
we introduce a soft-photon cutoff e on the minimum energy fraction of each 
photon, which is identical to the one used in the soft-photon contribution. 
We thus define: 

1 1 



Xi dxi / X2 dx2 5 ((1 — z) ~ (xi + X2)) 



(40) 



With the substitution xi {1 — z)-u and the abbreviation e = e/(l — 2), 
we have, after elimination of the trivial 5-function: 

l-e 

= {1 — z) du xi X2 ■ ■ ■ 

J L J xi={l-z)u,X2={l-z){l-u),... 



For the logarithmically (Lq) enhanced leading terms, we find the familar 
result: 

.2 



(2) 



log 



+ 



I 4- 

-4^^— lne + (1 + z) \nz - 2(1 - z) 
6(1 -z) + 

(2) 



^2 

-^0 



3+^^ 2 , .(1+^) 



1 - z 



In^ z + 4- 



1 - z 



In e 



2i±i!fln;2---21ne 
1-z V 2 



+ 



6(1 + 



S + z\ 2 . Al + zY, e 



1-z 



In^ z + 4- 



In- 



1-z 1-z 



with the leading-log radiator 
'1 + z^ 



3 \ 1 

21n(l-2;) - \xiz + - \ + -{1 + z)\nz-l + z 
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The analytical calculation of the remaining, non-enhanced terms is quite 
tedious, leading to lengthy expressions involving many dilogarithms and 
trilogarithms (see e.g., [P^]). 

As these terms also contain IR-divergent contributions, we shall pursue 
here the following approach. We analytically extract those terms in the inte- 
grand that either contribute to the infrared- divergence as e — or survive in 
the limit z ^ 1, before performing the integral over the remaining expression 
numerically. Besides, this separation improves the stability of the numerical 
integration. 

For the IR-divergent pieces of the non-enhanced terms we find 

= (^^\ = ^7du^ ^Ine (41) 

\(xiX2)V 1-^7 ''''^(l-^i) - 1-^ ^ ^ 

e 

We then decompose the infrared-finite pieces as follows: 

p(2),IR-fin./ N _ pz^l prem ( \ (Ar)\ 

nonlog V^/ nonlog ^ nonlogl^/ ' l^^y 

SO that -PnoSog(-^ = 1) = 0. The first term on the r.h.s., obtained in the limit 
z —>■ 1 reads: 



nonlog 



3XiX2 



Xi + 3X2 , , 3Xi + X2 , 

Inxi H ^ lna;2 



2 -^1 
1-Z)3 1-z 

ln(l - z) + 



xfxl " ' X1X2 



du- 



2m , 1 + 2m, , 

m M H — m(l — u) 



u{l — u) (1 — uY u 



= -y(3 + -^) , (43) 

where we neglected terms of order e. 

Finally, we plot in figure || the 'remainder' -PnoSogl-^)? which we evaluated 
numerically. 



B Non-factorizing quasi-collinear piece 

In the discussion of the contributions to the semi-coUinear piece, a non- 
factorizing part was split off in eq. (pSj). The function H appearing there is 
given by 

H{xi, X2) = Hi H H2-x{^ + ^ ] , (44) 

XiX2ri XiX2r2 \ri rt^ ) 
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prcm 
nonlog -30 




0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

Z 

Figure 2: The infrared-finite, 'normalized' part of the double-colhnear con- 
tribution, Pnonw(^)- 



with 



H3 



+ 



+ 



X 



— ^ cosip, 

iln^^ + ^+Li, (-^ 

2 X2X Q V 

2cos'?/' _ ,^ , 

Hi + — ^ln(2x2ri) 

1 + cos ip 

1 + 2 cos 



X 



1 



^ , , X2ri 

^ cosip, 

Xir2 

1 + cos ip 



1 + cos ip 
1 



1 + cos ip 
ln(?7 + X2ri + Xir2 cos ip) 



In 



1 + cos ip 

2co?,ib 

H2 + 



ln(?7 + Xir2 + X2ri cos ip) 



1 + cos ip 



ln[(l + cosip)x2ri] + 



1 + cos Ip 



In 



1 + cos Ip 



1 



1 + COS Ip 

1 + 2 cos tp 



1 + cos Ip 
I-X1-X2 , ri 2 



ln(r7 + X2ri + Xir2 cos ip) 
ln(?7 + Xir2 + X2ri cosip) 



1 -Xi 2 . 



The function S and expressions r] and cosip are defined in (|38|). For a deriva- 



tion and details we refer the reader to [BlI]. 
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